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Abstract 


A singularly perturbed time delay parabolic problem with two small pa- 
rameters is considered. The paper develops a finite difference scheme that 
is exponentially fitted on a uniform mesh in the spatial direction and uses 
the implicit-Euler method to discretize the time derivative in the temporal 
direction in order to obtain a better numerical approximation to the solu- 
tions of this class of problems. We establish the parameter-uniform error 
estimate and discuss the stability of the suggested approach. In order to 
demonstrate the improvement in terms of accuracy, numerical results are 
also shown to validate the theoretical conclusions and are contrasted with 
the current hybrid scheme. 
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1 Introduction 


We deal with the following class of singularly perturbed parabolic initial- 
boundary-value problems (IBVPs) on the domain D = Q, x (0,T], Q2 = (0,1): 


$-,yu(x,t) = ut — €Uze — pa(z,t)u, + (2, t)u = w(z, 6), 
u(a,t) = dp(2,t), (a,t) €T, = [0,1] x [-7,0], 
u(0,t) = drt), Ty = {(0,t):0<t<T}, 

J=@(t), Tr={(,t):0<t<T}, 


(1) 
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where w(x, t) = —c(x,t)u(x,t—7)+f (a, t), (z,t) € D. Here in this paper, T = 

T, Ul, UL, with the parameters ¢ and yz such that 0<¢ <1,0< <1, and 

T > Orepresents the delay parameter and the functions a(z, t), b(x, t), c(x, t), f(a, t), 
op(a,t), b(t), and ¢,(¢) are sufficiently smooth, bounded functions indepen- 

dent of ¢ and yu with 


a(z,t)>a>0, O(a,t)h>B8>0, c(a,t)>Vv>0, 
(x,t) € D = (0,T] x [0,1]. 


The type of singularly perturbed two-parameter problems changes depending 
on the values of the perturbation parameters € and py; for uw = 0, the problem 
is a reaction-diffusion problem, whereas, for pp = 1, it is a convection-diffusion 
problem. It is well known that due to the presence of layers, classical numeri- 
cal methods using a uniform mesh cannot properly approximate the exact so- 
lution when the parameter decreases unless a large number of mesh-intervals 
are utilized. However, even for lesser values of the perturbation parameters, 
one can overcome this difficulty by employing the fitted operator technique, 
which works without the prior location of the boundary layer. Time delay 
parabolic differential equations have recently attracted increasing amounts 
of attention due to their widespread use in many diverse application fields, 
including material science, biosciences, medicine, control theory, economics, 
and so on; see [20, 1, 10, 21, 23]. Many researchers have discussed the numer- 
ical results of the solutions of one-parameter singularly perturbed parabolic 
differential equations with time delay. For instance, one can refer to the 
articles by Das and Natesan [3], Gowrisankar and Natesan [6], Kumar [7], 
Woldaregay et al. [22], and Negero and Duressa [13, 14, 15, 16, 17]. 


In recent years, the development of a fitted numerical scheme for solving 
singularly perturbed time-delay parabolic problems having two parameters 
has received significant attention from a few authors. One such efficient fitted 
numerical scheme is an upwind difference scheme, which is proposed for solv- 
ing singularly perturbed time-delay parabolic problems having two parame- 
ters in [5] by Govindarao, Mohapatra, and Sahu. They constructed a method 
on the Shishkin type meshes (standard Shishkin mesh, Bakhvalov-Shishkin 
mesh) and proved that the method is first-order accurate. Negero [12] con- 
sidered the same problem in [5] and produced a second-order convergent 
scheme using an exponentially fitted cubic spline scheme. Prior to Negero’s 
[12] strategy, there were no developed numerical techniques for addressing 
two-parameter singularly perturbed time-delayed parabolic problems based 
on fitted operators. Kumar et al. [8] devised and analyzed a hybrid monotone 
finite difference scheme for singularly perturbed IBVPs of the form (1). In 
[8], a first-order uniformly convergent method is given for problem (1) using 
a hybrid monotone finite difference scheme on a rectangular mesh, which is 
a combination of a uniform mesh in time and a layer-adapted Shishkin mesh 
in space. There were no established numerical methods for dealing with 
two-parameter singularly perturbed time-delay parabolic problems based on 
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fitted operators prior to Negero’s [12] strategy. Thus, the main aim of the 
present study is to provide robust parameter uniform convergent numerical 
methods based on exponentially fitted for the solution of problem (1). 
Organization of the paper: In Section 2, the properties of the continuous 
solution are given. In Section 3, we describe the construction of an exponen- 
tially fitted finite difference discretization of problem (1). The stability and 
uniform convergence analysis of the suggested technique are given in Section 
4. Some numerical results that validate our theory are reported in Section 5. 
Lastly, in Section 6, we present the conclusion of the paper. 

Notations: The norm |]-|| is used to denote the maximum norm over the 
domain D, defined as ||g|| = maxp |g(zx,t)| for a function g defined on some 
domain D. In addition, C and its subscripts stand for positive constants 
independent of the perturbation parameters ¢, j1, and mesh sizes. 


2 Properties of the continuous solution 


The required compatibility conditions at the corner points are 


$v(0,0) = (0), 
{ du(1,0) = 64(0), (2) 
2e1(0) a ® pa(0, 0) 22010.) + b(0,0)¢4(0, 0) 


= —c(0,0)¢5(0,-—7) + f (0,0), 


2 
710) FO) — ua( 1,0) PE) 4 5¢4, 0)60(1,0) 


= —c(1,0)¢(0,—r) + f(1, 0), 


so that the data matches at the two corners (0,0) and (1,0). Let a, b, c, and 
f be continuous on a domain D. Then (1) has a unique solution u € C? (D) 
[9]. 


(3) 


Lemma 1 (Continuous maximum principle). Let ®(2,t) € C? (D)NC® (D) 
and ®(z,t) > 0 for all (x,t) ¢ =f,UT,UT,. Then $2.7 (x,t) > 0 in D 
gives ®(x,t) > 0, for all (x,t) € D. 


Proof. Assume (6*,¢*) € D such that ®(6*,¢*) = mini, yep (x,t) and 
@ (6*,¢*) < 0. Then, it is easy to verify that $.,,,®(0*,¢*) < 0, which is a 
contradiction. Thus, we have ®(a,t) >0 for all (x,t) € D. 


Lemma 2. The solution u(x,t) of the continuous problem (1) is bounded 
as 
|u (x,t) — dp (a, 0)| < Ct. 


Proof. Refer to [8]. 
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Lemma 3. The bound on the solution u(x,t) of the continuous problem (1) 
is given by 


ju(x,t))<C, (a, the D. 


Proof. Refer to [8]. 


Lemma 4 (Uniform stability estimate). Let u(x,t) be the solution of the 
continuous problem in (1). Then we have the bound 


||u(x, t)|] < B~* |lwl| + max {]dp| + max (|i , |r|) }- 


Proof. An application of Lemma 1 to the comparison function 


x* (a, t) =B~* |[gl| + max (|$o], (di + l@rl)) Eu(a,t), — (#, t) € D, 


yields the required estimate. 
Lemma 5. Let u(x,t) be the solution of problem (1), satisfying 0 <i+2j < 
4. Then u(a,t) satisfies the following bound: 

when ap? <n, 


s 2\J 
(4) (+) when ap? > en, 


Otsu 
| Ox? Ot) 


where 7 © ming t)¢ 


oi 
= 

8 

oe 
Nw 


Proof. Refer to [8]. 


3 Numerical scheme formulation 


3.1 Temporal discretization 


The time interval [0,7] is partitioned into a uniform step size as follows: 
QM = {tm = mAt,m=0,1,...,M,At=T/M}, T=ks, s=m,At, 


where k is a positive constant, m, is a positive integer, At is the time step 
size, and M is the number of mesh intervals. 

Hence, the problem (1) is discretized by using the implicit Euler method 
as follows: 
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um —Uu™ - 
“S @)_. (Use) ** (2) — pat" (2) 
Un" @) hi UP @) au), 
U™ (0) =diltm), OSMSM,xE Ms, 4 
U™ (1) = or (tm); O0<m<M,x €Qz, 
ies (x) = bb (2, tm+41) =(e-- 4) <ms< =i, rEs, 
where w™t1 (x) = —c™t (2)U™t1S (2) +. f™* (x), 0O< m < M,x € 9, and 


U™*1(x) is the semidiscrete approximation to the exact solution u(, tm41) 
of (1) at the (m+ 1)th time level. Then, let us rewrite (4) in the following 
operator form: 

Sh") = H (a, tm) 

um" (0) =i (tat), O<m<M, 

Ur (1) = abe (ln) 5 O0<m<M, cEM, 

um (a) = bp (2, tm4a)), —(s+1)<m<-l, rEM, 


where 
Sag aaa) =-—€ Cu (x)—pa™*? (x) ore (x)+q"t} (x) ymtt (x) 
and 


eta) _ role (x) _ goa (a) + f(g), 


1l<m<M, rEQ,, 


1 
for gg” (2) = a perl a). 


Lemma 6 (Semidiscrete maximum principle). Let y+! (x) € C?(D)N 
C® (D) . If y™*! (0) > 0, y™*? (1) > 0, and $@,y™"! (x) > 0 for alla e D, 
then y™+1(x) >0 for all z € D. 


Proof. One can prove this lemma by the same procedure as the proof of 
Lemma 1. 


Ou (x,t Z 
Lemma 7 (Local error estimate). Suppose ule) < C,(a,t) € Dx (0,7), 
0 <i< 2. The local truncation error defined as e41 = u(2,tm) —U™ (2), 


associated to scheme (5) satisfies 


llem+ill < C(At)?, m=1,2,...,M. 


Proof. See [2]. 


Lemma 8 (Global error estimate.). The global error E,,+1 is estimated as 
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|Em+ail] < C (At). 


Proof. See [3]. 


At the (n + 1)th time level, the characteristics equation of the homoge- 
neous part of the differential equation (5) can be 


ed? (x) + pat (x) A(x) — (om (x) + x) = 0. (6) 


Then, the roots of (5) are 


2€ 
_ geth(a = qmtl r 2 * 
ete ee i( ue) a 


1 
where o* = b’*! (x)+—. From these roots, it is possible to see the boundary 


layer behavior of the solution in the neighborhood of « = 0 and x = 1. Let 
00 = —MaXze/0,1] A1 (7) and @1 = minz¢{9,1] A2 (7). Then we have two cases 


i) When EU, ase 0, oo © 01 = y/ 2, where 0 < of < 0". 

ii) When ie + 0,as p> 0, = £05 and @; = 0, where 0 < 93 < parte), 
Next, we give the semidiscrete bound of the solution U™*1(2) of the problems 
in (6). 


Lemma 9. [8]For a fixed number 0 < p < 1 and for a certain order k, the 
solution U™() of (5) satisfies the following derivative bound 


d'U™ (cr) 


dx* 


eC (1 + gp te Pet 4 grees) ) , for 0<i<k. 


3.2 Fully discrete problem 


In this section, we fully discretize the problem under consideration via an 
exponentially fitted finite difference scheme for space derivative discretiza- 
tion. On the space domain [0,1], we introduce the equidistant meshes with 
uniform mesh length h such that 


ON = {ay = nh, n=1,2,...,N,29 =0,2n =1,h =1/N}, 


where hf is the step size, and N is the number of mesh points in the space 
direction. Using the theory applied in the asymptotic method developed 
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in [18], we develop an exponentially fitted numerical scheme to solve the 
singularly perturbed BVPs in (6). In the considered case, the boundary layer 
is on the left side of the domain, so for the singularly perturbed problem of 
(6), the zero-order approximation asymptotic solution is given as 


U1 (a) UP (0) + (G1 (tm) — UN (0)) exp { A (x 2) ac} 


€ 
+ O(e), 
(7) 


where Uj”*'(x) is the solution of the reduced problem in (6) obtained by 
setting ¢ = 0 written as 


m d m m m m 
pal" (a) Uy (x) — g(x) Ug" (2) = G(x), 


UPS OO) = oi Ger) 


(8) 


Taking Taylor’s series expansion for a(a,t) about x = 0 and taking their 
first terms, (7) gives 


U2) UB) + (6 (maa) = UB*O)) ex { = (HO?) a} 
+O(e). 


At the mesh x,, = nh, (9) becomes 


U™+ (nh) =Ug"t* (nh) + (1 (tm4i) — UZ"t*(0)) exp {- (2) inn) 
+O(e). 
Therefore, 


lim U4 (nh) = Ug"**(0) + (01 (tm) — U3"**(0)) exp {na (0) np}, 


(10) 
where p = Bh 
Now, we consider the derivative approximation of the problem in (1) and 


(2) as 


U, — Un- Un —U, n ~~ Un— 
DU, = = 4 D*U,, = pel Aa D°U, = a ne =. and 
Uns — 2Un + Un 
47— _ n+1 n n-1 
P= 2 
and 
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eo (p,€,2) D}D-U™*? (an) + pa™** (a) DOU™™ (an) 


- (11) 
—q me (i) U Per (tn) =G (2s begat) ’ 
where o (p,€, 4) is a fitting factor. 
Multiplying (11) by A and evaluating the limit as h > 0 give 
| (PEs) f pr m m 
[202 (a 20") 
. (12) 


1 
+ 5a" (nh) (Oey! — UR") = 0. 


Substituting (10) into (12) and taking a(z,t) = a constant with some manip- 
ulation give the fitting factor as 


git 
o(p,€,4) =a! (0) $ coth (2). 


For the variable fitting factor, we define as 


pa™ 


On (p,£,H) = a"™*# (am) £ coth ( (13) 


Hence, using (12) ,the resulting finite difference scheme can be given as 


1 
via Chea = (= (p, E, 7) par) Shae 


ue h2 2 
(=: (P, E, 7) a) ue 
h (14) 
EOn (p,€, 4) 1 m+1 m+1 
( h2 F glen Uns 
=m 
subject to the following conditions: 
UR adi Gaels 0<m<M, 
Unt! = ¢r(tmti); O<m<M, (15) 
Uiagtass) = bp (Rr teh) te S Oo —(9+1) <m< —l, 
where 


Hy" = H (2n,tm+1) 
1 


At U™ (tn) + eae Ga prrnie). 


The schemes in (14) and (15) can be rewritten as 


MS eee a ee Gee ae, (ie) 


n 
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where 
Em+ — Eon tens) = ; gmt, 
prt = 2600 (EH) 4 cr’, 
gmt _ £On ers LL) a Suan, 


1 = m 
Be = ag (t,) +e" (2, )U™- PO (@,) — faa). 


From the entries E7’t!, F*1!,G™*1 of tridiagonal system of (16), it is evi- 
dent that E™*1 < 0,G™t! <Oand E™14 Fpmt1 + Gmt! s 0. Thus the 
system is an M-matrix, and therefore its inverse exists, and it is positive. 
Hence, the tridiagonal system in (16) can be easily solved by any existing 
methods. 


4 Stability and uniform convergence analysis 


Lemma 10 (Discrete maximum principle). Assume that 7™*1 is any mesh 
function that satisfies yj't' > 0, wt! > 0, and that $N,” is the discrete 
operator of (16). Then $¥Mym+t > 0, for 1 <n < N —1, implies that 
wm > 0, forO<n<QN. 


Proof. Refer to [12]. 


Lemma 11 (Uniform stability estimate for discrete problem). Let U?"t! be 
any mesh function such that Uj"*! = 0,Un't! =0 on 0<n< N. Then 
max |§) @Urt| 


jog] < Pee 


+ C max {|di(tm41)|5|@r(tm4i)[} 5 


1 
where gt! = me +O tite) 2g > 0; 


Proof. Refer to [12]. 


Theorem 1. Let U (2n,tm+41) be the continuous solution of (1) and (2) and 
let U"*! be the approximate solution of (16). Then, for sufficiently large N, 
the following error bound holds: 


| ae (U (an, tm-+1) — ur) <—ON 


Proof. Consider the error bound in the spatial direction as 
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sai (U (tn, tm41) — iia | 


= (8° U Gating oF 


E (ee ee (—) + pa (2a) oy (2a) 


— {20 (p,¢,u) DtD-URY + par Durty} 


(17) 


d d 
ca (p,6,4) ($5 - DYD-) Up + want! (5 — p°) op 


< 
dx? 


< 


mt+1 
e Ge (2n) = coth (PE“ , (®n) | Dp ye 


~ + 


ae + 41 
—— DTD )U™ 


part (= _ °) ym 
xv 


Now, (17) becomes 


[sa (UT (@nstmaa) — Un") | 


277m+1 Arrm+1 3 TTM+1 
f a Cen? + Cpe, 


dx? 


< Cuh? 


Using Lemma 9, we have 


GU Gained) UP) 


<C ph? (1 fs a mda al oh wy 2e-watl-2) | 


1 Ch? € (1 +4 ee lad de cage ae) 


+ (1 ae ee ee fs ajtenae-a) : 


As € > 0 both wy ’e7’1*™ and wy 'e~’2—-*m) approach zero for 0 < i < 4. 
Thus, we obtain the following error bound: 


|" (U (anstm4i) — Ut)| < CN, 


since h = N7-!. 


Under the hypothesis of Lemmas 11 and 10, the following error estimate 
holds: 


t —U™+|) << Ch=CN-?. 18 
pmax, |U (tn; m+1) — UR |<c C (18) 
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Theorem 2. Let u(z,t) be the exact solution of (1) and (2) and let U™*? 
be the numerical solution of (16). For the discrete scheme, there exist a 
constant C independent of ¢,h and At such that 


O<n<Ndem<M | 24: (ara treet) _ ue <C (N~? + (At)) 3 


for sufficiently large N. 


Proof. The result follows from the error estimate given in Lemma 8 and 
Theorem 1. 


5 Numerical results 


In this section, we illustrate the proposed scheme using two numerical ex- 
amples of the form given in (1). We investigate the theoretical results in 
this paper by performing experiments using the proposed scheme. The exact 
solution of these two examples is not known. Thus, we use the double mesh 


principle to evaluate maximum absolute errors EX ue and the corresponding 
N.M ag 


order of convergence pz; 


EN.M 
N,M _ m+1 Im+1 NM _ ef 
Eo OP Ua |, Bey = 10g, (aoe 


ey max | 
0<n<N,0<m<M 
SoS Eb 
From these values, we obtain the e-uniform error EN” and the corresponding 
e-uniform order of convergence p™™ as 


N/M _ N,M N/M _ Not 
EN = Maxo<n<No<m<mu ES’ and p'” = logs (fea : 


where U”*? is the numerical solutions obtained by using N x M mesh inter- 
vals in space and time direction with mesh size h and At, respectively. 


Example 1. Consider the problem 


Ou Oru Ou 7 9 
a a (1 + t) + u(x,t) = —u(2,t —7)+162* (1-2), 
(x,t) € (0,1) x (0, 2], 
with 


Example 2. Consider the problem 
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Ou OPu 2, Ou 
De epg Ht a(l—2) +h) So + (1+ Sat) ue, t) 


=-u(z,t—T)+2(1—2) (e'-1), (x,t) € (0,1) x (0, 2], 


with 


Table 1: Maximum pointwise errors (EN, 7) and rate of convergence (ps i) for Example 


1. 


p= 10-4 N=32 N=64 N=128 N=256 N=512 
et M=16 M=32 M=64 M=128 M=256 
10-9 5.751 6e-03 3.6382e-03 2.1286e-03 1.1627e-03 6.0947e-04 
0.66074 0.77332 0.87243 0.93185 = 
10-2 1.0422e-02 5.449 1e-03 2.7875e-03 1.4101e-03 7.0919e-04 
0.93554 0.96705 0.98317 0.99155 7 
10-4 1.0658¢-02 5,531 2e-03 2.8189e-03 1.4231e-03 7.1502e-04 
0.94627 0.97246 0.98610 0.99298 = 
10-6 1.0663¢-02 5.5328e-03 2.8193e-03 1.4232e-03 7.1504e-04 
0.94653 0.97267 0.98620 0.99304 - 
10-8 1.0664e-02 5.5339e-03 2.8202e-03 1.4237e-03 7.1533e-04 
0.94638 0.97250 0.98615 0.99296 = 
10-19 1.0664e-02 5.5339e-03 2.8202e-03 1.4237e-03 7.1533e-04 
0.94638 0.97250 0.98615 0.99296 7 
10-¥ 1.0664e-02 5.5339e-03 2.8202e-03 1.4237e-03 7.1533e-04 
0.94638 0.97250 0.98615 0.99296 - 
ENM 1.0664e-02 5.5339e-03 2.8202e-03 1.4237e-03 7.1533e-04 
pNM 0.94638 0.97250 0.98615 0.99296 . 


2 


° 


Numerical Solution - U 
2 


Numerical Solution - U 


2 
PY 


an 


(a) (b) 


Figure 1: Surface plot of the numerical solution for Example 2 with N = 256, M = 128, 
ae=10-!,4=10-!7 be = 10717, = 10-1 
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N,M 


Table 2: Maximum pointwise errors (Ez), 
1. 


) and rate of convergence (pe ap for Example 


y= 10-' Number of mesh intervals N=M 
et 32 64 128 256 512 
10~° 3.6218¢-03 2.1253e-03 1.1619e-03 6.0925e-04 3.1225e-04 
0.57919 0.87118 0.93138 0.96433 - 
10-2 5.4110e-03 2.7780¢-03 1.4077e-03 7.0860¢-04 3.5549e-04 
0.96185 0.98071 0.99030 0.99516 : 
10-4 5.5307e-03 2.8187¢-03 1.4231e-03 7.1501¢-04 3.5838¢-04 
0.97243 0.98599 0.99300 0.99647 - 
10-° 5.5315¢-03 2.8189¢-03 1.4231e-03 7.1502¢e-04 3.5839e-04 
0.97254 0.98610 0.99298 0.99645 - 
10-8 5.5315¢e-03 2.8189¢-03 1.4231e-03 7.1502¢e-04 3.5838e-04 
0.97254 0.98610 0.99298 0.99645 : 
10-1° 5.5315¢-03 2.8189¢-03 1.4231e-03 7.1502¢e-04 3.5838¢-04 
0.97254 0.98610 0.99298 0.99645 - 
10-12 5.5315¢-03 2.8189¢-03 1.4231e-03 7.1502¢e-04 3.5838¢-04 
0.97254 0.98610 0.99298 0.99645 : 
ENM 5.5315e-03 2.8189e-03 1.4231e-03 7.1502e-04  3.5839e-04 
pNiM 0.97254 0.98610 0.99298 0.99645 - 


Table 3: Maximum pointwise errors (EN, Ae) and rate of convergence (ps iM) for Example 
>: 

=10-7 N=32 N=64 N=128 N=256 N=512 

el M=16 M=32 M=64 M=128 M=256 

10-9 2.1475e-04 1.091 2e-04 5.4962e-05 2.7578e-05 .3813e-05 
0.97674 0.98941 0.99492 0.99749 - 

10-? 2.1465e-03 1.1561e-03 6.0053e-04 3.0592e-04 .5440e-04 
0.89272 0.94496 0.97308 0.98648 - 

10-4 2.6764e-03 1.4488¢-03 7.5345e-04 3.8424e-04 1.9401e-04 
0.88544 0.94327 0.97150 0.98588 - 

10-6 2.6771e-03 449 1e-03 7.5407e-04 3.8484e-04 .9445e-04 
0.88551 0.94239 0.97044 0.98486 - 

10-8 2.6771e-03 .4490e-03 7.535 1e-04 3.8423e-04 .9401e-04 
0.88561 0.94336 0.97166 0. 98584 F 

10-19 2.677 1e-03 A490e-03 7.5351e-04 3.8423e-04 .9401e-04 
0.88561 0.94336 0.97166 0. 98584 - 

1o-P 2.6771¢-03 4490e-03 7.5351e-04 3.8423¢-04 .9401¢-04 
0.88561 0.94336 0.97166 0. 98584 - 

ENM 2.6771e-03 1.4491e-03 7.5407e-04 3.8424e-04 1.9445e-04 

pNM 0.88551 0.94239 0.97269 0.98261 - 
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‘i ) and rate of convergence (ps Fras ) for Example 


j= 10-'? Number of mesh intervals N=M 
et 32 64 128 256 512 
or? 1.3372e-04 6.1251e-05 2.9166e-05 1.4212e-05 7.0123¢-06 
1.1264 1.0704 1.0372 1.0191 - 
0-? 1.1701e-03 6.0326¢-04 3.0645e-04 1.5447e-04 7.75606¢-05 
0.95578 0.97713 0.98833 0.99394 - 
o-4 1.4466¢e-03 7.5262e-04 3.8382¢-04 1.9386¢-04 9.7408¢-05 
0.94267 0.97149 0.98541 0.99290 - 
0-6 1.4522e-03 7.5525e-04 3.8509¢e-04 1.9444¢e-04 9.7702e-05 
0.94321 0.97176 0.98587 0.99287 - 
0-8 1.4522e-03 7.5527¢e-04 3.8510¢e-04 1.9445¢-04 9.7705e-05 
0.94318 0.97176 0.98583 0.99289 - 
g-10 1.4523e-03 7.5527¢e-04 3.8510¢e-04 1.9445¢-04 9.7705e-05 
0.94328 0.97176 0.98583 0.99289 - 
o-P 1.4523e-03 7.5527e-04 3.8510¢e-04 1.9445¢-04 9.7705e-05 
0.94328 0.97176 0.98583 0.99289 - 
ENM 1.4523e-03 7.5527e-04 3.8510e-04 1.9445e-04 9.7705e-05 
pNM 0.94328 0.97176 0.98583 0.99289 - 
Table 5: Comparison of uniform error (EN’™) for Example 1. 
p=10% N=32 N=64 N=128 N=256 N=512 
et M=8 M=16 M=32 M=64 M=128 
Proposed method 
10-4 1.9859e-02 1.0660e-02 5.5318e-03 2.8190e-03 1.4232e-03 
10-® 1.9905e-02 1.0684e-02 5.5439e-03 2.8245e-03 1.4251e-03 
10-8 1.9905e-02 1.0684e-02 5.5440e-03 2.8252¢€-03 1.4262e-03 
10-19 1.9905e-02 1.0684e-02 5.5440e-03 2.8252e-03 1.4262¢-03 
1o-!? 1.9905e-02 1.0684e-02 5.5440e-03 2.8252e-03 1.4262¢-03 
Method in [8] 
10-4 4.3705e-2 1.6704e-2 7.3802e-3 3.7406¢-3 1.8967e-3 
10-® 4.3471e-2 1.6596e-2 7.3290e-3 3.7218¢-3 1.8873e-3 
10-8 4.3429e-2 1.6573¢-2 7.3303e-3 3.7211le-3 1.8870e-3 
10-10 4.4343e-2 1.6572e-2 7.3303e-3 3.7211le-3 1.8870e-3 
10-12 4.4343e2 1.6572e-2 7.3303e-3 3.7211e-3 1.8870e-3 
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Table 6: Comparison of uniform error (E™’™) for Example 2. 


p= 10-* N=32 N=64 N=128 N=256 N=512 
et M=8 M=16 M=32 M=64 M=128 
Proposed method 
o-4 4.5627e-03 2.6876¢e-03 1.4564e-03 7.5785e-04 3.8653e-04 
0° 4.5254e-03 2.6603e-03 1.4402e-03 7.4904e-04 3.8203e-04 
o-8 4.5254e-03 2.6603e-03 1.4402e-03 7.4904e-04 3.8195e-04 
Oras 4.5254e-03 2.6603e-03 1.4402e-03 7.4904e-04 3.8195e-04 
o-? 4.5254e-03 2.6603e-03 1.4402e-03 7.4904e-04 3.8195e-04 
Method in [8] 
o-4 1.1161e-2 5.1087e-3 2.4749e-3 1.2214e-3 6.0706e-4 
o-° 1.1008e-2 5.0450e-3 2.4437e-3 1.2073e-3 6.0036e-4 
o-8 1.0941e-2 5.0426¢e-3 2.4442¢e-3 1.2071e-3 6.0016e-4 
g-t0 1.0940e-2 5.0428¢e-3 2.4442¢e-3 1.2071e-3 6.0016e-4 
o-? 1.0940e-2 5.0428¢e-3 2.4442e-3 1.2071e-3 6.0016e-4 


Table 7: Comparison of uniform error (EN) for Example 1. 


p= 10-9 N=32 N=64 N=128 N=256 N=512 

ei) M=8 M=16 M=32 M=64 M=128 
Proposed method 

10-4 1.9853e-02 1.0658¢-02 5.5313e-03 2.8189e-03 1.4231¢-03 
10-6 1.9856¢-02 1.0659e-02 5.5315e-03 2.8189e-03 1.4231e-03 
10-8 1.9856¢-02 1.0659e-02 5.5315e-03 2.8189e-03 1.4231e-03 
10-16 1.9856¢-02 1.0659e-02 5.5315e-03 2.8189e-03 1.4231e-03 
10-12 1.9856¢-02 1.0659¢e-02 5.5315e-03 2.8189e-03 1.4231e-03 
Method in [8] 

10-4 4,.3708e-2 1.6705e-2 7.3807e-3 3.7407e-3 1.8967¢-3 
10-6 4.3816e-2 1.6749e-2 7.401 7e-3 3.7489e-3 1.9008¢-3 
10-8 4,.3817e-2 1.6750e-2 7.4019e-3 3.7490e-3 1.9008¢-3 
10-10 4.3817e-2 1.6750e-2 7.4019¢e-3 3.7490e-3 1.9008e-3 
107? 4.3817e-2 1.6750¢-2 7.4019¢e-3 3.7490e-3 1.9008e-3 


Table 8: Comparison of uniform error (EN) for Example 2. 


p= 10-9 N=32 N=64 N=128 N=256 N=512 

el M=8 M=16 M=32 M=64 M=128 
Proposed method 

10-4 4.5499e-03 2.6744e-03 1.4477e-03 7.5307e-04 3.8398¢e-04 
10-6 4.5651e-03 2.6830e-03 1.4523e-03 7.5527e-04 3.8513e-04 
10-8 4.5652e-03 2.6831e-03 1.4523e-03 7.5529e-04 3.85 14e-04 
10-10 4.5652e-03 2.6831e-03 1.4523e-03 7.5529e-04 3.8514e-04 
107? 4.5652e-03 2.6831e-03 1.4523e-03 7.5529e-04 3.85 14e-04 
Method in [8] 

10-4 1.1053e-2 5.0755e-3 2.4578e-3 1.2132e-3 6.0309e-4 
10-® 1.1046e-2 5.0765e-3 2.4625e-3 1.2161e-3 6.0456¢e-4 
10-8 1.1100e-2 5.0838¢-3 2.4627e-3 1.2161e-3 6.0457e-4 
10-16 1.1093e-2 5.0782e-3 2.4639e-3 1.2162e-3 6.0457e-4 
1o-? 1.1092¢-2 5.0775e-3 2.4640e-3 1.2162e-3 6.0457e-4 
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Max. Abs. Error 
Max. Abs. Error 


Figure 2: Log-Log plot of the maximum error on left (a) for Example 1 with » = 10-4 
and on right (b) for Example 2 with = 107+. 


We have demonstrated maximum pointwise errors ee, M) and the rate 
of convergence (p: oe M) for Example 1 using scheme (16) by fixing = 1074 i 
Table 1 and p = 10-2? ia Table’? with various-yalues of. Similar ly, Tables 3 
and 4 have presented the result obtained for Example 2. The results given in 
Tables 1-4 clearly indicate that the proposed numerical method is accurate 
of order O (N 24 At), which approves the hypothetical result predicted in 
the theory. Numerical solutions obtained by the presented numerical scheme 
(16) for Example 2 have been shown in Figure 1 (a), (b), and it shows the 
effects of the two parameters ¢ and yu on the steepness of the layer of the 
solutions. From Figure 1 (a), we confirm the nonoccurrence of both left and 
right boundary layers near x = 0 and x = 1 for yp + 0 as € becomes large. 
Similarly, from Figure 1 (b), we confirm the occurrence of left boundary 
layers near x = 0 for uw = 1 as € becomes small. The graphs between N and 
maximum pointwise errors of Examples 1 and 2 are plotted as the log-log 
scale, respectively, in Figure 2 (a) and (b). From these two graphs, one can 
observe that the numerical scheme converges ¢-uniformly as the perturbation 
parameter goes very small. The comparison of our numerical results with 
that of [8] is presented in Tables 5-8. From these tables, we can confirm the 
improved accuracy of our proposed numerical method. 


6 Conclusion 


A singularly perturbed parabolic differential equation exhibiting boundary 
layers was considered. The considered problem contains two small pertur- 
bation parameters multiplied by the highest order derivative a term of the 
equation and a large delay parameter on the time variable. An exponentially 
fitted operator numerical scheme was proposed for solving the problem. First, 
the equation was approximated by equivalent singularly perturbed parabolic 
partial differential equations using the implicit Euler method in the time 
direction. Inducing an exponential fitting factor for a term with the per- 
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turbation parameter ¢ and determining its value, a fully discrete numerical 
scheme was developed using implicit Euler in temporal discretization and 
the central finite difference method for spatial discretization. The uniform 
stability and uniform convergence of the scheme were established. It was 
shown that the scheme is accurate and converges uniformly with the order 
of convergence O (N~? + (At)). 


Acknowledgements 


The author is grateful to his anonymous referees and editor for their con- 
structive comments. 


References 


[1] 


wy 


Asl, F.M., and Ulsoy, A.G. Analysis of a system of linear delay differ- 
ential equations, J. Dyn. Sys., Meas., Control, 125 (2) (2003), 215-223. 


Clavero, C., Jorge, J. and Lisbona, F. A uniformly convergent scheme 
on a nonuniform mesh for convection-diffusion parabolic problems. J. 


Comput. Appl. Math. 154(2) (2003), 415-429. 


Das, A. and Natesan, S. Uniformly convergent hybrid numerical scheme 
for singularly perturbed delay parabolic convection-diffusion problems on 
Shishkin mesh, Appl. Math. Comput. 271 (2015), 168-186. 


Epstein, LR. Delay effects and differential delay equations in chemical 
kinetics, Int. Rev. Phys. Chem. 11 (1) (1992), 135-160. 


Govindarao, L., Sahu, S.R. and Mohapatra, J. Uniformly convergent 
numerical method for singularly perturbed time delay parabolic problem 
with two small parameters, Iran. J. Sci. Technol. Trans. A Sci. 43(5) 
(2019), 2373-2383. 


Gowrisankar, S. and Natesan, S. ¢- uniformly convergent numerical 
scheme for singularly perturbed delay parabolic partial differential equa- 
tions, Int. J. Comput. Math. 94 ( 2017), 902-921. 


Kumar, D. A parameter-uniform scheme for the parabolic singularly 
perturbed problem with a delay in time, Numer. Methods Partial Differ. 
Equ. 37 (1) (2021), 626-642. 


Kumar, S. and Kumar, M. A robust numerical method for a two- 


parameter singularly perturbed time delay parabolic problem, Comput. 
Appl. Math. 39(3) (2020), 1-25. 


Iran. J. Numer. Anal. Optim., Vol. 13, No. 4, 2023,pp 627-645 


Negero 644 


[9] Ladyzhenskaia, O.A., Solonnikov, V.A. and Ural’tseva, N.N. Linear 
and quasilinear equations of parabolic type.(Russian) Translated from the 
Russian by S. Smith Translations of Mathematical Monographs, Vol. 23 
American Mathematical Society, Providence, R.I. 1968. 


[10] McCartin, B.J.  Discretization of the semiconductor device equa- 
tions,” New problems and new solutions for device and process mod- 
eling”, Boole, (1985), 72-82. 


[11 


Miller, J., O’Riordan, E., Shishkin, G. and Shishkina, L. Fitted mesh 
methods for problems with parabolic boundary layers, Math. Proc. R. Ir. 
Acad. 98A (1998), no. 2, 173-190. 


[12 


Negero, N.T. A uniformly convergent numerical scheme for two param- 
eters singularly perturbed parabolic convection-diffusion problems with a 
large temporal lag, Results Appl. Math. 16 (2022), Paper No. 100338, 15 
Pp. 


[13] Negero, N.T. and Duressa, G.F. A method of line with improved accuracy 
for singularly perturbed parabolic convection-diffusion problems with large 
temporal lag, Results Appl. Math. 11 (2021), 100174, 13 pp. 


[14 


Negero, N.T. and Duressa, G.F. An efficient numerical approach for 
singularly perturbed parabolic convection-diffusion problems with large 
time-lag, J. Math. Model. 10(2) (2022), 173-110. 


[15] Negero, N.T. and Duressa, G.F. Uniform convergent solution of singu- 
larly perturbed parabolic differential equations with general temporal-lag, 
Iran. J. Sci. Technol. Trans. A Sci. 46(2) (2022), 507-524. 


[16] Negero, N.T. and Duressa, G.F. An exponentially fitted spline method 
for singularly perturbed parabolic convection-diffusion problems with large 
time delay, Tamkang J. Math. (2022). 


[17] Negero, N.T. and Duressa, G.F. Parameter-uniform robust scheme for 
singularly perturbed parabolic convection-diffusion problems with large 
time-lag, Comput. Methods Differ. Equ. 10 (4) (2022), 954-968. 


[18 


O’Malley Jr, R.E. Introduction to singular perturbations, Applied Math- 
ematics and Mechanics, Vol. 14. Academic Press [Harcourt Brace Jo- 
vanovich, Publishers], New York-London, 1974. 


[19] Roos, H.G. and Uzelac, Z. The SDFEM for a convection-diffusion prob- 
lem with two small parameters, Dedicated to John J. H. Miller on the oc- 
casion of his 65th birthday. Comput. Methods Appl. Math. 3(3) (2003), 
443-458. 


[20 


Tikhonov, A.N. and Samarskii, A.A. Equations of mathematical physics, 
Courier Corporation, 2013. 


Iran. J. Numer. Anal. Optim., Vol. 13, No. 4, 2023,pp 627-645 


645 


A robust uniformly convergent scheme for two parameters... 


[21] Van Harten, A. and Schumacher, J. On a class of partial functional dif- 
ferential equations arising in feed-back control theory, Differential equa- 
tions and applications (Proc. Third Scheveningen Conf., Scheveningen, 
(1977), pp. 161-179, North-Holland Math. Stud., 31, North-Holland, 
Amsterdam-New York, 1978. 


[22] Woldaregay, M. M. Aniley, W. T. and Duressa,G.F. Novel numerical 
scheme for singularly perturbed time delay convection-diffusion equation 
Adv. Math. Phys. (2021), Art. ID 6641236, 13 pp. 


[23] Wu, J. Theory and applications of partial functional differential equa- 
tions, New York, Springer, 119, (2012). 


How to cite this article 

Negero, N.T., A robust uniformly convergent scheme for two  pa- 
rameters singularly perturbed parabolic problems with time delay. 
Tran. J. Numer. Anal. Optim., 2023; 13(4): 627-645. 
https://doi.org/10.22067 /ijnao.2023.80721.1214 


Iran. J. Numer. Anal. Optim., Vol. 13, No. 4, 2023,pp 627-645 


